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Abstract 


A technique for measuring smoke plume of large indu%trlal 
sources observed by satellite using LARSYS is proposed. A 
gausslan plume model Is described. Integrated in the vertical, 
and Inverted to yield a form for the lateral diffusion coeffi- 
cient, Ky. Given u, wind speed; y^, the horizontal distance 
of a line of constant brightness from the plume symmetry axis 
a distance downstream from reference point at x=x 2 . y“0, 
then Ky * uy ^/2xj^ln(x2/x2^) . The technique Is applied to a 
plume from a power plant at Chester, Virginia, imaged August 
31, 1973 by ERTS I. The plume bends slightly to the left 4.3 
km from the source and estimates yield Ky of 28m2s“ near the 
source, and 19 m^s”^ beyond the bend. Maximum groif ' concen- 
trations are estimated between 32 and 64 ugm"^. Foisting met- 
eorological data would not explain such concentration . 
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List of Symbols 1 

» distance down wind from plume effective source 

■ horizontal distance from plume symmeti ’ axis | 

■ vertical distance from plume symmetry axis = 

“ mass per unit time of effluent emitted 

“ mean scalar wind speed 

“ Concentration of effluent as mass per unit volume 
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- lateral effective diffusion coefficient 
(l^t“^) 

■ vertical effective diffusion coefficient 

(l2t”l) 5 

■ (x^ + + 2^)'^ 

■ columnar loading in vertical in mass per unit 
area 

» number of i*-^ pixel in scan line j on ERTS 
image 

■ number of scan line on ERTS image 

■ distance in meters from p^j to pj^^^ in config- 

uration space 

- compass bearing from ground point Pj^]^ to 
in configuration space at latitude = 37° 


Introduction 

Remote sensing techniques have been applied to the study 
of smoke plumes from fixed Industrial point sources In an 
effort to verify numerical smoke plume dispersion models * In 
previous papers^’ ^ the value of such studies was discussed, 
together with techniques that would apply to aerial photo- 
graphs of plumes. In brief, such studies can aid in the dev- 
elopment and refinement of plume dispersion modelr that com- 
bine stack and effluent parameters, with meteorological, or 
climatological data as input, and that yield ground level de- 
position patterns for use in air pollution control schemes 
and in environmental impact studies of new industrial sources. 
In addition, if such models exist, they can be inverted so 
that comparison with Images of smoke plumes can yield valuable 
insights to regional pollution transport, wind trajectories, 
and vertical stability at the time the image was acquired. 

This is obviously usefu 1 during pollution episodes, or during 
field experiments. 

The curi^ent study extends previous work by changing 
from aircraft to satellite platform (ERTS-I, recently renamed 
lANDSAT I) Imagery, and by applying computer techniques to the 
analysis of the digiti 2 ed Imagery itself. To do this, use was 
made of the LARSYS software developed by the Laboratory for 
Applications of Remote Sensing (LARS) at Purdue University, 
Lafayette, Indiana, and of a remote terminal located at the 
National Aeronautics and Space Administration's Langley Re- 
search Center (LRC) , Hampton, Virginia. The terminal gave 
direct access to the LARS computer and databank at Purdue, and 
was jointly administered by LRC and by Old Dominion University 
of Norfolk, Virginia. I i 
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For the current study a detailed analysis was made of 
the Image of a smoke plume (ERTS E- 1404-15190) as It appeared 
(Cf. Fg 1) about 10:00 EST, August 31, 1973. The plume 
arises from a steam-electric generating plant In Chester, 

I Virginia, operated by the Virginia Electric Power Company 

(VEPCO), and It can be followed on the image toward the north- 
west about 18 kilometers. An estimate was made of the lateral 
diffusion coefficients as measured from ^he Image. Some 
anomalous meteorological conditions were examined. Including 
an interesting large scale bend In the plume, and the fact 
that the wind field as reported by the nearest available met- 
eorological stations was at odds with the observed behaviour 
of the plume. Finally, the plume model was run using the data 
acquired fvom the image, and estimates of ground deposition 
were made. 


i 



Fig. 1 Smoke plume at Chester, Virginia, 10:00 EST, as seen 
by ERTS MSS 4 (LARS channel 1), 31 August 1973. North is to- 
ward the top. The light vertical region branching from the 
plume middle is a ground feature (Army Munitions Depot). 
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Analytical Techniques 

The plume model used will be briefly described outlining 
its required input parameters, its limiting assumptions and 
the forms of output it can produce. Next the type of data 
available through LARS will be described. Finally, two pro- 
blems in relating the model to the LARS satellite data, spa- 
tial rectification of the LARS printout and computation of 
lateral diffusion coefficients, will be outlined. 


Plume model . The computer plume model attempts to solve 
a turbulent diffusion equation for a time Invariant point 
effluent source in a wind field. The input data required is 
in two categories — stack and meteorological parameters. 

Stack parameters, available from the State Air Pollution Con- 
trol Board, Include stack height and diameter, and the velo- 
city, temperature, and mass flux of the effluent. Meteorolog- 
ical parameters incluLe mean wind direction and speed, surface 
temperature and lapse rate. Standard empirical formulas can 
be applied to those data to determine an "effective stack 
height" for the plume. ^ The plume’s effective source point is 
usually somewhat above the true stac3c height due both to the 
vertical effluent momentum and to its buoyancy as it enters 
colder ambient air. 


The mathematics of the turbulent diffusion equation, and 
its approximate solution, sometlm e s called the gausslan plume 
model, is covered in a variety of texts.^>^ Some of its sal- 
ient features will be described here. First, it assumes that 
the situation is steady state, both for stack and meteorolog- 
ical parameters. This assumption is frequently valid, but ' 
often variability in winds causes the solution to be valid 
only on a time average basis. Second, it is assumed that in 
any vertical plane normal to the mean wind vector the concen- ■ 
tration distribution of effluent is elliptically distributed • 
about the plume’s symmetry axis with the major and minor axes 
aligned in horizontal and vertical directions. Third, it is 
assumed that the concentration falls off in any such plane 
along any vertical or horizontal line as a gaussian distri- 
bution centered at the major or minor axis of the ellipse. 

This assumption is only valid if a set of effective turbulent 
diffusion coefficients can be determined Such that the lateral 
coefficient, Ky, and the vertical coefficient, K^, are both 
independent of spatial and temporal coordinates and such that | 
the downwind coefficient, is zero to first order. Bound- j 
ary conditions are applied forcing concentrations to zero at , 
large distances from the source. No et fluent sinks are con- j 
sidered. The mass flux out the stack, 0 equals ’Fu | 
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integrated over any plane normal to the plume's symmetry axis. 
Concentrations then can be described as 

T(x,y,z) = (Q/AirrCKyK^)^) exp ((-n/4x) (y^/Ky + z^/K^)) 

While Kg, can be estimated from meteorological conditions 
(lapse rate, wind sheer, etc.) the determination of Ky Is al- 
most Impossible from a priori considerations. 

Several output forms are available from the model. One 
plots the concentration along any plane orthogonal to any co- 
ordinate. A second provides an output that is called the col- 
umnar loading, the definite integral of f from minus infi- 
nity to plus infinity over z, which is related to brightness 
observed from a satellite. This is given by; 

<I'(x,y) S (Q/(-\iruKyx)^) exp (~uy2/4xKy) 

Finally, the model is capable of superimposing multiple plumes 
of differing stack parameters and locations, and of modifying 
downwind meteorological conditions , 

LARS data . While the LARS system has a variety of soft- 
ware, only certain ones of the less sophisticated algorithms 
were employed. These included the LINEGRAPH, COLUMNGRAPH, 
and PICTUREPRINT programs. The LINEGRAPH program will take 
any scan line of any spectral band, and will prepare a graph 
of the radiance as a function of pixel (picture element) num- 
ber along the scan line. The COLUMNGRAPH feature plots rad- 
iance versus scan line number for a given pixal number. Fin- 
ally, PICTUREPRINT prepares a greyscale map of a region in a 
single ERTS spectral band with an option to assign up to 16 
symbols to 16 different grey scale domains specified by the 
user. When used in conjunction with auxiliary programs (HIS- 
TOGRAM' and GRAPHISTOGRAM) it is possible to density slice the 
map in ,such a way that background is completely suppressed 
(i.e. printed as a blank), while a feature of interest (e.g. 
a smoke plume) is printed in several symbols representing dif- 
fering regions of brightness or reflectivity. Figure 1 shows 
such a printout though photoreduction for publication has 
made it hard to distinguish the different symbols. Lines of 
constant brightness can then be projected on the horizontal 
plane, and the first two graphing programs can be used to det- 
ermine inflection points on the spatial brightness curve, 
yielding lines of maximum contrast that can be superimposed 
by hand on the grey scale map. 
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Fig. 2 LARSYS PICTUREPRINT of smoke plume. 

ERTS MSS 4 (LARS 1) 

LARSYS PICTUREPRINT Spatial Rectification . The grey scale 
maps printed by the LARS System (LARSYS) are spatially dis- 
torted. ERTS Is In a near polar orbit at 99° to equator. 

Thus, the ground trucK Is not along lines of constant longi- 
tude. The MSS scans perpendicularly to the ground track of 
the satellite. Thus the MSS Image Is a parallogram projected 
on the earth's surface. In any given line a pixel of a cer- 
tain number represents a point slightly west of the point due 
south of the similarly numbered pixel In the preceding line. 
The computer's printer Inherently prints pixels of Identical 
number directly beneath one another. Moreover, the number of 
typed symbols per line, and lines per Inch do not correspond 
directly to the actual ground scene spacing of areas repre- 
sented by pixels and lines. To correct these problems, three 
simple formulas were developed so that true distances and 
compass bearings could be determined from a LARSYS PICTURE- 
PRINT. Select two pixels Plj and Pkl on scan lines J and 1, 
then the actual distance In meters, between the two 

points on the ground Is given by 

^Ijkl " (79m)((Xi-Xj)^ (Pij-p^l)2(0.494))‘5 
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The azimuthal angle, in degrees measured the east 

of 'north between the two pixels is then 

^ijkl ” tan“^((p^j-Pjj.j) (0.70)/(X2“^j)) 

Finally, we will be interested in measuring horizontal plume 
widths on grey scale maps along lines perpendicular to the 
downwind symmetry axis of the plusae. Thus, given two pixels, 
P£j and pj^j^, on that symmetry Sxis, we would like to find 
another pixel, in any line, 1^^ of our choosing that re- 

presehts a point on a line passing through p. . perpendicular 
in configuration space to the pluftie. axis. This is given by 

Pm - Plj.+ (2.02)((X„-Ap(Xj-Xj)/(Pij-p^p) 

Theory for determination of the lateral diffusion co- 
efficient . In a previous publication^ measurements of lateral 
diffusion coefficient by eye from aerial photos were reported^ 
The eye detects the "shape" of the plume as the line of maxi- 
mum contrast between plume and background, and not as a line 
of constant brightness. In that case it was reported that 
the shape of the line of maximum concentration is parabolic, 
with the apex at the source, and in accordance with the model 
it satisfies the formula: 

y = (2Kyx/u)^ . 

. Through the LARSYS PICTUREPRINT program It is possible to 
determine lines of constant brightness. We may assume these 
are related, to lines of constant columnar loading, as long as 
the plume is not opaque, and is of sufficiently low concen- • 
tration for multiple light scattering events to be insigni- 
ficant. The shape of such a line can be derived from Eq2 by 
setting $(x2,0)=#(x2^,y2) . The result is 

M}(TruKyX2)“^ 2 %Q(7ruKyX3^)“^xp(-uy|/4xxKy) • . 

This can be solved for the contour of lines of constant ^ 

brightness. 

yX ((2xxKy/u) In (x 2 /xx))^ . | 

Such lines form closed contours (Cf Fg. 3) and if x's and y’s/ 
can be determined experimentally from measurements on a gray | 
scale map one can solve for Ky ' 

O 9 ’ 

Ky“nyx/(2xx In (x^/xx)) j 
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Thus, In theory, the lateral diffu- 
sion coefficient can be measured. 


Practical considerations . Due 
to the fact that a real, turbulent 
pluine is being interpreted with the 
aid of an Idealized time averaged 
model, it is insufficient to use one 
set of data points (x’s and y's) to 
determine the lateral diffusion 
coefficient. A graplilc overlay 
technique is used to smooth data sc 
that a series of data points are 
used to determine an average value 
of Ky. In essence, a plume s 3 Tnmetry 
axis is drawn on the gray scale map 
by eye. Second^ a series of lines 
perpendicular to it are marked out 
using the spatial rectification Fg 3. Lines of constant 
procedure. Then a smooth out- brightness as seen from 

line of the plume’s lines of above, 

constant brightness are sketched. 

Widths are measured as a function 

of downwind distance, Ky's are computed and averaged. These 
are then used as input to the plume model. 

One other point is that different ERTS spectral bands 
give different results. This is to be expected since the 
ability of smoke particles to scatter light is strongly de- 
pendent on particle size and optical wave length. ERTS band 
7 (LARS. band 4) in the near infra-red (0.8 - l.lura) cuts 
through haze readily, and frequently a plume that is very 
visible in another band is almost undetected in band 7. 

This may '.prove to be a useful tool for obtaining a crude 
size distribution for plume aerosols. 

Analysis of the Chester, Virginia, Plume 

A beautiful example of an industrial smoke plume was 
imaged by ERTS-1 on Friday, August 31, 1973, arising from 
an oil fired power generating station at Chester, Virginia. 
All input information needed to run the model could be 
obtained, and the image was taken with ideal atmospheric 
conditions. Thus, it was decided to analyze the case in 
detail. The source of the plume will be described, as will 
the meteorological conditions. A peculiar bend in the 
plume will be briefly discussed, then lateral diffusion 
coefficients will be estimated. Finally, the results jof • 
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the plume model run using tb'^e values of Ky will be dis- 
cussed^ 

Plume source . The ,;ower plant at Chester, Virginia is 
located on the south bank of the James River, 17 km NNE of 
Petersburg, and 19 km SSE of Richmond. The plant has six 
smoke stacks. The height, h; internal diameter, I.D.; 
effluent efflux velocity, v; .sj ’luent temperature, T; and 
rate of fuel (No. 6 oil) consu^iiption, f, are tabulated 
below. 


Stack 

h(m) 

I.D.(m) 

v(ms"’^) 

T(Oc) 

f (kgs"^) 

•1 

61 

3.7 

12 

197 

5.91 

2 

61 

3.7 

15 

174 

7.96 

3 

6l 

4.0 

14 

142 

9 . 75 

, 4 

61 

4.0 

22 

162 

15.0 

5 

61 

5.2 

23 

130 

29.3 

6 

61 

6.1 

33 

128 

58.9 


Meteorological conditions . The analysis of meteorolo- 
gical conditions at the time the plume was imaged presents 
a conundrum, because the reports of the nearest stations 
indicate a situation at odds with the observed plume behavi- 
or. Eastern Virginia was under the influence of a weak 
high pressure system with associated light breezes. Winds 
reported at Richmond were calm at 9:58 EST and reached 7 
knots from 200° (SSW) at 10:56 EST. The satellite imaged 
the scene at 10:14 EST (15:14 Z) . From sunrise until 9:58 
EST, the temperature had risen rapidly from 69°F to 89°F at 
a rate of about 5°F/hour, but from 9:58 EST until 10:56 EST 
the temperature rose only 2°F, and during the following 
hour it failed to rise at all. Such behavior would indi- 
cate that a weak diurnal Inversion persisted over the 
region until the mbming insolation caused sufficient heat- 
ing of the earth's surface to warm the lower layer of air. 
When the inversion broke up, sometime between 9:58 EST and 
10:56 EST a much thicker mass of air had to be heated 
resulting in the decreased rate of temperature increase. 
Though no radiosondes are flown from Richmond, the neatest 
ones (from Wallops Island, Virginia, and from Washington, 

D. C.) do indicate the presence of such an inversion that 
morning . ' 

1 

The presence of this inversion, and its breaking up, is 
in accord with the observed plume behavior, but the wind 
directions are not. All surrounding stations reported winds 
out of generally southwesterly directions, yet the plume 
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clearly was heading off toward the northwest, thus being 
blown by southeasterly winds. The significance of this fact 
will be discussed In the conclusions, while the accord with 
the Inversion will be discussed In, the next section. 

Bend In Plume , An obvious feature of the plume is that 
it bends. Starting from the smirce, the plume heads 4.3 km 
at an azimuth of 305°, then turns westward at 272°. At 
first this bend was assumed to indicate the rise of the 
plume to a higher altitude with differing flow direction, 
but this peculiar bend to the left was not consistent with 
the. tendency of winds at higher altitudes to blow at an 
angle to the right of the surface winds. The idea was pro- 
posed to the authors by Herman Wobus and Earl Kindle that 
the bend did not represent a spatial variation in the wind- 
field, but rather a temporal one. Prior to approximately 
one half hour before the image was acquired, the wind did 
blow from the east carrying the plume away from the source 
at 272®. Then it shifted to a more southeasterly wind that 
carried the effluent subsequently discharged off in bear- 
ing of 305°. This wind shift is associated with downward 
momentTom transfer after the breaking up of the inversion. 

But significantly, the new wind direction carried the older 
plume, still with its symmetry axis oriented at 272® bodily;, 
off at 305°. This changed the Ky value obtained for por- 
tions of the plume farther from the source than 4.3 km 
from the value obtained closer to the source. Thus, two 
different lateral diffusion coefficients were to be ex- 
pected for the two regimes. Wobus predicted, in fact, that 
a smaller Ky in the more distant portion of the plume could 
allow higher concentrations to occur at the bend point (a 
virtual source) than would be experienced closer to the 
real stack, Superficially, an examination of the greyscale 
maps indicated a region of enhanced brightness was indeed 
present, but it has so far been impossible to definitively 
establish that it is not related to a ground feature, a 
large military depot, that is unfortunately right beneath 
the bend. A LARSYS CLASSIFYPOINTS program did remove the 
background feature leaving the bright downstream region, 
but its photometric interpretation is difficult. When the 
two values of Ky determined from the imagery were inserted 
in the plume model, it did display the feature Wobus pre- 
dicted. » 

Determination of lateral Diffusion Coefficients . The 
coefficients determined by the method above can be deter- 
mined as a function of wind speed. Between the source and 
the bend care had to be exercised to not allow the multiple 
stacks to confuse the issue by causing an apparent 
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broadening of the plume that would be Interpreted as an ex- 
cessively large value of Ky. Thus, only a limited number of 
points of adequate downstream distance could be used. The 
mean value was determined to be Ky=ll,6 (ra) u i 0.6 (m) u. 

In the region beyond the bend more numerous measurements lead 
to values of Ky = 8,0 (m)ut 1.6 (m)u. 


Wind speeds were now estimated by noting that the bend 
in the plume is 4.3km from the source. The meteot^logical 
data Indicated the inversion was breaking up around the time 
of the last observation at 9:58 EST. Assuming the wind 
shifted about a half hour before the satellite pass, the wind 
speed would have to be approximately 8.6km per hour, or to 
one significant figure 6 knots. This makes the values of dif- 
fusion coefficient between source and bend Ky * 28 m^s 
This value is probably overestimated due to apparent broaden- 
ing of multiple sources. Assuming a similar windspeed beyond 
the bend, the value would bw Ky = 19 m2s“l. Both these values 
arr?. of course based on crude assumptions. 


Model predicted groun d concentra tlons . Inserting the 
above values into the plume model, together with stack and 
meteorologidal parameters, the following output (Fg 4) was 
produced. Concentrations at ground level are printed on a 
binary exponential scale, such that a 2 indicates 2^ iiom-3 

ugm”3. 


ugm 


and a 5 would indicate 2- 
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Maximum ground concentration Inside tho region marked 5 would 
be greater than 5'^ or 32 ugm“^, but less than 2® or 64 ugm”3. 

Conclusions 

.• The application of satellite imagry to the analysis of < 
smoke plumes has shown that useful data can be acquired that 
could not be assembled in any other way. First, anomalously- 
high particulate pollution levels that would potentially be ■ 
recorded by surface air quality stations would have been I 
attributed to the wrong cause on the test day simply because; 
the nearest meteorological stations reported different wind | 
directions. Consequently, local air pollution models for ; 
urban regions that depend on a limited number of reporting { 
stations should be checked against remotely sensed data to | 
Identify the extent and and frequency of such anomalous 
behavior. Second, the smoke plume itself records a temporal, 
history of the day’s meteorological conditions, filling in 1 
gaps in locally recorded data. Finally, if emission parame- j 
ters are known satellite data can lead to reasonable esti- 
mates of surface loading due to a plume simply by seeing its ‘ 
apparent geometric shape. 
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